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ABSTRACT 

A generalized v/aterjet propulsion system in a sub- 
cavitating hydrofoil craft is considered for design and 
performance estimation. Independent system parameters of 
jet velocity ratio, inlet velocity ratio, and nacelle inlet 
diameter to maximum diameter ratio are varied to search for 
the optimum system. The optimum system is defined as the 
minimum total propulsion system weight. The optimization 
scheme utilizes a directed search without the calculation 
of derivatives, and was chosen for its simplicity, versa- 
tility, and rapidity. The ducting system is divided into 
components for head loss calculations. For a given design, 
total head losses in the duct are computed by means of 
experimental data and empirical equations, enabling pump, 
reduction gear, and prime mover design to be completed. 
Results for a sample craft are included which indicate that 
hydrofoils should be designed about the gas turbine, due to 
the discrete power levels that are available. Suboptimizing 
the nacelle design may lead to an overall less optimum 
system design. A general lack of information was noted on 
nacelle, cascade corner, and three dimensional diffuser 
design and performance prediction. A FORTRAN computer 
listing and flow charts are included. 
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CHAPTER 1 - INTRODUCTION 

The use of waterjets for the propulsion of marine 
vehicles at conventional ship speeds has not been consid- 
ered comparable to that of propellers. Although the con- 
cept is quite old, the disadvantages of waterjets have not 
been overcome and the advantages not realized. At the 
speeds considered, the low system efficiency, large system 
weight and relatively sophisticated technology required 
were sufficient to virtually exclude waterjets from com- 
petition with propellers. With the advent of high speed 
craft however, this means of propulsion is receiving a 
renewed interest as a practicable system. Among the pos- 
sible advantages that may be realized are (refs. 1,2): 

a. Reduced noise and vibration 

b. Fewer system components 

c. Elimination of underwater appendages 

d. Elimination of complex transmission machinery 

e. Limited draft applications 

f. Simplification of steering devices 

For hydrofoil craft, the waterjet looks particularly 
attractive. It neatly avoids the problem of conventional 
propellers in power transmission (refs. 3,4). For this 
reason, waterjet propulsion development to-date has been 
primarily concerned with hydrofoil applications. In 



addition, new concepts in ships, such as surface effect 
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ships and captured air bubble craft, are now including con- 
sideration of waterjets as a possible propulsion system. 
With this surge of interest, the need to go beyond the 
practical "Will it work?" engineering design method has 
become paramount. The question has become, "What is the 
best waterjet propulsion system for the craft under consid- 
eration?". Answering this question requires a strictly 
formulated methodology previously not required in make-it- 
work technological levels. No longer is it satisfactory 
to find some point, any point, above the minimum acceptable 
level. The scope has been enlarged to pinpointing the 
optimum within the entire domain of feasibility. The 
immensity of this task coupled with its relative lack of 
need has hampered the implementation of a methodology. 
Hydrofoils have demonstrated the need, while computers 
provide the invaluable aid for the required computations; 
therefore, an optimization methodology has been developed. 
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CHAPTER 2 - DESIGN METHODOLOGY 



To initialize the design method# a clear definition 
of the desired objective must be given. In the present 
instance# the definition is the answer to the question# 
"\’7hat is meant by an optimal water jet propulsion system?". 
The choice of an optimum must adequately describe the 
system without bias. It is readily apparent that whatever 
optimum is chosen# it will have a strong effect on the 
methodology employed. The current most commonly used 
descriptor is the propulsive efficiency which is defined 
as the ratio of the power required to the pov/er output of 
the plant (ref. 5). Efficiencies suffer from differences 
in definition and thus vary in their relative importance to 
system optimization. Minimum power required to propel the 
craft at the desired speeds is some.times utilized as an 
optimization objective with the advantage of being readily 
calculable and understood. A disadvantage is that a gas 
turbine# which is a commonly used prime mover on account of 
its low weight to power ratio# operates most efficiently at 
its maximum power level. Hence# there exists a natural 
trade-off between the most efficient propulsion and the 
margin required for growth in the ensuing years of the 
craft's life. Minimum power calculations are certainly 
necessary but provide a somewhat questionable optimim. In 

cin attempt to avoid vagaries of definition and usage and 
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to remove doubts of sufficiency, the goal of optimization 
or the system "cost" may be defined differently. The 
"cost" of any system should reflect the disadvantages of 
that system v;ithout weighting one disadvantage over 
another. Two disadvantages of waterjet propulsion systems 
mentioned previously are a large system weight, in part 
due to the entrained water, and a low overall system 
efficiency. Inefficiencies are, however, necessarily 
manifested in increased weight, be it added required 
fuel, another pump stage, or larger ducting. Recall, too, 
that hydrofoils are weight limited vehicles. A unique 
number that reflects these realizations is a suitable 
definition of the weights attributable to the propulsion 
system (ref. 5). The propulsion system weight shall be 
defined as the sum of all weights of the structural com- 
ponents, prime movers, fuel, pumps, gearbox, ducting, and 
entrained water above the mean waterline. If, for a given 
craft, the displacement is considered constant, then minimi- 
zation of the propulsion system weight is equivalent to 
maximizing the craft payload. Thus, a propulsion system 
shall be defined as optimal if, in addition to fulfilling the 
propulsion requirements, it adds the minimum weight of 
any propulsion system to the vehicle. 

The generalized waterjet propulsion system for 

hydrofoils utilized in this study is depicted in Figure 1. 

It should be noted that this configuration permits 
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TOP VIEW 




FIGURE 1 GENERAL DUCTING CONFIGURATION 



considerable variation by the changing of a few lengths 
and angles of the ducting. The water enters the nacelle 
under ram pressure and is ducted up the strut and into the 
pump where the energy input is converted to kinetic 
energy in the jet. Thus the thrust reaction is generated, 
providing the craft's propulsion. In equation form, 

T = pQ(Vj cos (B) - Vq) (1) 

In any optimization procedure, the formulation of the 
method requires the knov/ledge of the parameters to be 
varied. The number and selection of these parameters 
contributes significantly to the applicability of the 
method. The parameters must be necessary and sufficient 
system descriptors, and it is desirable to utilize the 
minimxam number required. If M operations are performed in 
the evaluation of one design point and N independent para- 
meters describe the system, then the total number of 
operations required increases as M^. The sufficiency of the 
description of the system is clearly a choice of adequacy. 

The variables that most strongly affect the resulting system 
should be chosen as system parameters while other variables 
that produce only minor variations over a wide range of 
values on the essential design should be defined, or even 
somewhat arbitrarily chosen, to suit the needs of the design. 

Examination of equation (1) reveals that for a given 

craft and associated mission, only the flow rate and jet 

velocity are vinknown. The one -dimensional continuity 
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equation for an incompressible fluid 

= A 2 V 2 = Q 

indicates that choice of jet velocity will determine the 
jet area. The first non-dimensionalized system parameter, 
then, shall be Vj/Vo, defined at the cruise point. Knowing 
the flow rate from the jet, the same reasoning may be applied 
to the inlet which results in the second system parameter of 
Vj^/Vq, where Vj^ is defined in the cruise condition at the 
point of minimum area in the nacelle. The third and 
last system parameter, for reasons that will be explained 
in Chapter 3, shall be Dj^/Dj^, where Di is the diameter of 
the inlet at the point of . 

Finally, the minimum amount of information needed to 
specify the craft and its mission must be determined. . A 
typical drag to craft speed curve is shown in Figure 2. 

The "hxmp" at the lower speeds represents the take-off 
point, i.e., when the hull lifts free of the water and the 
craft is foilborne, and generally occurs at speeds that are 
approximately 50% of the cruise speed. The sizing of the 
propulsion system then is dependent on both the cruise and 
take-off conditions. The minimum information required is 
the following: 

a. Displacement 

b . Range 

c. Prime mover 

d. Craft Sneed, at take-off and cruise 
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FIGURE 2 TYPICAL HYDROFOIL LIFT TO DRAG VS. SPEED CURVE 



e. Craft Drag, at take-off and cruise 

f. Depth of submergence of foil 

g. Distance of strut from transom 

h. Distance of pump exit from transom 

i. Height of pump centerline above keel 

Although the prime mover is assumed to be a gas tur- 
bine, knowledge of the particular type of gas turbine is 
important in consideration of both fuel and power re- 
quirements. It is readily seen that the last four input 
requirements define the approximate configuration while 
the first five items describe the craft size and mission. 

Knov;ing the required ingredients and the desired re- 
sult leaves only the recipe to be set forth. For the 
traverse of an ideal fluid between two points in a closed 
conduit, Bernoulli's equation may be applied such that 

Pi + + pgzj^ = P2 

where 2 indicates a station downstream of station 1. 

For a real fluid with losses this becomes 

Pi + ^ ?2 " ^^Pi2 

where ~ looses between stations 1 and 2. 

Although this equation can be applied only along a stream- 
line and hence cannot be used across the pump, it is 
applicable to the ducting on either side of the pump. The 
continuity equation provides the information that the 
ducting velocities are determined by geometric configura- 
tion. Starting from a known point, e.g. the water surface, 
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and following the water through the duct leaves only the 
losses as unknown. Estimation of the losses may be accom- 
plished based on duct size and shape (ref. 8) . This yields 
the needed information for the pump design. 

Following the flow diagram of Figure 3, and recalling 
the previous discussion, the choice of the system parameters 
by the technique described in Appendix A determines the in- 
let and outlet geometry. Knowledge of the jet area determines 
the flow rates at both cruise and take-off. Water is in- 
gested into the nacelle at velocity V^. The parameter 
sets the amount of diffusion possible prior to entry into the 
strut elbov;. Determination of the nacelle length provides 
both a drag estimate and the head loss in the nacelle duct 
as described in Appendix B. The strut elbow will exhibit 
large losses due to the low aspect ratio (width/depth) of 
the duct as it conforms to the shape required for the strut 
cind due to the relatively large velocities. The use of 
guide vanes will minimize the losses by increasing the 
aspect ratio. Assuming the inlet and outlet areas are the 
same, the losses in the elbow may be calculated by the method 
described in Appendix C. The strut diffuser area ratio is 
set to avoid cavitation in the dov/nstream ducting and to 
reduce the inflow velocity to the pump. Sizing the internal 
flow area permits an estimate of the exterior dimensions 

which are sized up to cavitation limits. 

-19- 



For more information. 



/ INOUT t 



I 

I 

I ^ 

Specify ci-AfT 
CONf IGLRATICN *NO 
NISSI IN 



--0 I , n 4» — >1 
P4RMS I 



I SfT SYSTIH I 

I PARA^eiPfS I 



ARE VALJES Of * 
SYSTcN PA<<AfETCRS • 
WITHIN RTUNOS? • 



SEE NJT£ 
AfiOVE 



I CALCULATE JET ANO | 
I INLET AREAS I 



I 

OETERNINE FLOW 
RATE AT CRUISE 
ANO TAKE-DFF 
CONDITIONS 



I 

NACELLE I 07 

I SIZE NACELLE, I 

1 OETERNINE DRAG I 

I AND MEAO LOSS I 



WILL THE 
NACELLE 
CAVITATE? 



OETERNINE HEAD 
LOSS IN STRUT 
ELOJW 



I SIZE STRUT, I 

I oetfrnine drag I 

I AND HEAD nSS I 



NO •!$ PREVIOUS* 

• ORAf, estimate 

•ACCEPTABLE?* 



I 



12 



I oetcrhine oucting I 

1 CONf IGURATtON ANO 1 
I HEAD LOSS UP TO I 
I PUMP I 

• — 

1 



I SIZE NOZZLE ANO I 
I DETERNINE HEAD | 
I LOSS I 



I oeterhine size I 

I AND TYPE Of PUMP I 



HAS OUCTING 
CONFIGURATION 
LSED ACTUAL PUMP 
length? 



I 



SEE NOTE 
ABOVE 



GEAR 

CALCULATE GEAR 
RATIO, GEARBOX 
WEIGHT, ANO 
REQUIRED SHP 

1 



•IS THE gear* 
RATIO LESS 
• THAN ONE? • 



YES 

• -♦ 

1 

I 

I 






• IS THE • NO 
REQUIRED SHP *-♦ 
•AVAILABLE? • I 



. 02 . 

PARMS 



SUM ALL I 

PROPULSION SYSTEM | 
WEIGHTS I 



IS THE PRESENT 
DESIGN BETTER 
THAN THE PREVIOUS 
BEST? 



SEE NOTE 
ABOVE 



I 

02 



I SAVE parameters I 
I ANO VALUE t 

I 

I 

I NOTE 2A 
*•••••••*•• 

• HAS THE MINIMUM • 

• WEIGHT SYSTEM • 

• BEEN ATTAINEO? • 



• NO 

SEE NOTE 

ABOVE • 1 

• • I 

• • I 



IVES 



I 

02 



I WITH OPTIMUM 
I SYSTEM 

I PARAMETERS. 

I RECONSTRUCT 

I OPTIMUM SYSTEM 

I 

I 

OUTPUT I 2 

/ 

z z 

Z OUTPUT SYSTEM Z 
Z OATA Z 

Z Z 



I 2B 

• halt • 



FIGURE 3 DESIGN METHODOLOGY FLOW DIAGRAf4 

- 20 - 



see Appendix D. The drags attributable to the strut and spray 
are estimated. The newly calculated nacelle, strut, and 
spray drags are compared to the original estimates. If the 
deviation is excessive, the drag estimate is revised and 
the problem restarted. Otherwise, the head losses in the 
strut due to friction and diffusion are calculated. . Entry 
into the hull from the strut is made via a 90® elbow for 
which the losses are calculated as before. It should be 
noted that this elbow will have a more favorable aspect 
ratio as a result of the strut shape and will require fewer 
guide vanes. Coupled with the decreased velocity, the 
losses will show a decrease from those due to the strut 
elbow. Determination of the ducting required, as described 
in Appendix E, permits estimation of the ducting losses 
therein. Summation of all the losses prior to the pump and 
knowledge of the inflow velocity provide the cavitation 
information in the form of net positive suction head (NPSH) 
to the pump. Adding the nozzle losses permits calculation of 
the punp head required. The method for determining tlie 
nozzle losses may be found in Appendix G. This is sufficient 
to determine the size and type of pump needed, the details 
of which may be found in Appendix F» The specifications of 
the prime mover are then utilized to specify the gearbox 
parameters of gearbox ratio and power which in turn indi- 
cates the SFC required in the cruise condition. Since the 
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amount of fuel required for the endurance specified is 
calculable, the entire propulsion system is fully defined 
in dimensions, configuration, weights, and particular com- 
ponents. Summing all the component weights yields the 
desired system weight. This concludes the evaluation of 
a propulsion system weight for the chosen system parameters. 
The optimization routine utilizes the weight information, 
selects a new set of system parameter values and the 
cycle begins anew. Ultimately, the minimum weight system 
is located and the process completer. 
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CHAPTER 3 - DISCUSSION OF RESULTS 



A sample output is depicted in Figure 6. The first 
page is predominantly a repetition of the input data and 
is printed out upon call to the program. The second and 
third pages are the resulting optimized system data that is 
printed for each configuration of pumps and gas turbines 
considered. Hence the designer can see at a glance the 
pertinent information required for a preliminary design as 
well as compare the alternative configurations that might 
be feasible and competitive designs. 

There are a number of features within the results that 
are not immediately obvious. First and foremost is the 
problem of results exceeding a system limitation. For 
example, having required cruise shaft horsepower exceeding 
the maximum normal horsepov^rer of the engine. Within the 
program during the optimization, such designs are, of 
course, rejected. For the print out, it was considered that 
knowledge of the resulting "optimum" was desirable, how- 
ever erroneous it may be. One of the forcing influences on 
this decision was the variability of the total drag estima- 
tion and thus the flow rate. If the optimum system re- 
sulted in a value that was at a boundary, a slight change 
in, say, flow rate might cause the value to exceed the 
boondary and hence result in a rejected design. In other 
words the optimum system is not identically reproducible. 

The significance of this fact will be amplified in the 
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conclusions section. This lack of reproducibility results 
primarily from the iterations for drag estimation. Each 
new estimation is compared to the previous estimation and 
accepted if within 5% of the former value. Otherwise the 
flow rate is adjusted for the new drag and the nacelle and 
strut are redesigned. In the case of the 400 ton craft 
used, 5% of the total drag is on the order of 4000 pounds 
force in the cruise condition. Thus, starting from the 
same system parameter values but v;ith different initial 
drag estimates may result in a total drag differing by 5% 
in the output v^7ith corresponding differences in the system. 

Secondly, the resulting design is obviously constrained 
by the ducting configuration used. Hence direct comparison 
with existing craft is not valid until the configuration 
is modified to conform to that of the craft. 

At the inception of this project, it was believed 

that the duct area ratio of the strut was a system parameter 

to be varied independently. Clearly a lower bound on this 

parameter is that amount of diffusion required to avoid 

cavitation in the hull ducting. It was felt, however, 

that the pump, gear and fuel v>?eight calculations might 

reflect a lower overall weight if the pump inflow rate was 

reduced. The results did not justify this assumption. The 

system unerringly chose the strut with the minimum required 

diffusion, which meant minimum. weight. The principal 

reason for this result is that the pump design is in terms 
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of total head. Hence the pump effectively never "sees" the 
amount of static head presented to it. The pump is affected 
only by the differences in ducting losses which are re- 
flected in both the NPSH and pump head. Thus the net effect 
is of secondary importance and is not significant enough 
to change the area ratio of the strut. For this reason the 
area ratio of the strut ducting was rejected as a system 
parameter and is presently determined as required by 
cavitation limitations or pump inlet area. 

The initial nacelle design method chose the largest 
inlet to maximum diameter ratio permissible by exterior 
cavitation. This ratio was tempered by the amount of 
diffusion required to avoid cavitation in the strut elbow. 
Assuming the same inlet velocity ratio, would more diffusion 
than the minimvim required permit an -overall lesser weight 
system? In other words, does the optimxm system require 
the optimum nacelle or might an off design nacelle be 
used resulting in a lesser weight system? The use of the 
nacelle diameter ratio, DIDM, as a system parameter pro- 
vides the opportunity to investigate this question. The 
current results are somewhat inconclusive due to the influence 
of other factors. Tentatively, it is suggested that for some 
designs, the use of DIDM as a system parameter is justi- 
fied and results in a better system design. It is 
recommended that investigation of this question be pursued 
for determination of what constitutes the basis of nacelle 
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design in regards to system performance. 

Perhaps the most significant conclusion is that the 
choice of engine is extremely important in the resulting 
optimum system. The first reason for this is that the 
gas turbine is presently available only at discrete power 
levels. Since fuel weight is large relative to the other 
system v;eights, the operation of a gas turbine at a point 
significantly off its design point will result in a large 
increase in fuel weight and thus system weight. The net 
effect is that, given .the opportunity to avoid the penalty 
in fuel weight, the system tends towards the maximum SHP 
permitted by the engine, and will drive the jet velocity 
ratio as high as possible to achieve this SHP. This also 
results in a low flow rate, low drag and low propulsive 
coefficient. If the engine is slightly large for the craft, 
then the results v;ill show a high jet velocity ratio, SHP 
per engine at the maximum permissible level, which may be 
limited by take-off or cruise SHP, and a low propulsive 
coefficient. If the engine is so large that there is no 
power limitation, the jet velocity ratio will be very 
high and the SHP per engine relatively low resulting in a 
very low propulsive coefficient due to large total SHP 
and a- high system weight. The effect of an engine that is 
barely adequate is a lowered jet velocity ratio in an effort 
to reduce the SHP requirements. The trade-off here is the 



large flow rate required to compensate for the low jet 
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velocity ratio and the large drag that results from 
attempting to accomodate the flow rate in the nacelle and 
strut. Thus SHP is being forced up and the system may never 
achieve an acceptable horsepower level. 

. This power level sensitivity is so strong that it is 
abundantly clear that hydrofoils should be designed 
around a particular gas turbine. Choice of a desired 
displacement and the use of an adequate prime mover may 
result in a poor design and inefficient use of the gas 
turbine. Choice of the gas turbine and permitting the 
displacement to be chosen at the point of best engine 
operation v;ill result in a more efficient and more effective 
design. Thus it is recommended that another version of this 
program be created that chooses the best displacement 
given the engine. It is also suggested that perhaps the 
concept of hydrofoil craft design be revised to accomodate 
the fact that the engine choice determines the best design. 
Figure 7 shows the minimum weight ratio for one LM2500 
and one pump with a lift to drag ratio of 11.3. The 
results show the optimum weight ratios with and without 
fuel do not coincide. The propulsive coefficient is drawn 
to show the trade-off between maximizing SHP and maximizing 
the propulsive coefficient. 

The choice of an engine also affects the reproducibility 
of the optimum system design. If the engine is smaller 

than desirable, VJVO is driven as high as the SHP permits. 
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When the optimum system parameters have been found and the 
optimum system is to be regenerated, the design process is 
reinstituted. During this final design, iterations in 
craft drag are again required to determine the final flow 
rate and size of the system. Thus the net result may be 
a system slightly different from the one originally found 
with least weight. And in the case of a system already on 
a boundary, the resulting system may exceed the limits 
placed on the optimization process. It is to be expected 
that the deviation will be small but nonetheless present. 

The results in such a case should be interpreted as optimal 
at the boundary and not an unsatisfactory design. If the 
deviation is significant auid the other system values 
reflect it, the interpretation should be that the system 
was unable to be optimized within the limits. The most 
common example of this occurs when the cruise or take-off 
SHP of the engine is insufficient to meet the craft re- 
quirements. 

The effect of range on the system design is shown in 

Figure 4. The LM 2500 is a power limited engine for the 

four hundred ton craft at forty-two knots. The FT4C-2, 

however, is so large that maximum power is never attained in 

this craft. The effect of lower SFC and operation at 

maximum power favors the LM 2500 at the longer ranges. As 

it should be, the weight ratio without fuel is essentially 

constant for both engines, although the weight ratio for the 
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LM 2500 is higher, presumeibly on account of the power limita- 
tion. Similarly, the jet velocity ratio remains constant 
as does the flow rate, but the inlet velocity ratio tends 
to increase with increasing range. The nacelle diameter 
ratio, on the other hand, shows no consistent pattern. 

The effect of variation of displacement by 10% is shown in 
Figure 5 for the FT4C-2 for one and two pumps. 

The results show a surprising sensitivity to the 
starting values of the system parameters. When arbitrarily 
deviated by twenty per cent from the present values, the 
results were unpredictable. The current values, however, 
have shov/n good stability with test craft ranging from 
fifty to seven hundred fifty tons and from forty to fifty 
knots. The reason for this apparent instability is un- 
known. The step sizes appear to be acceptable at any size 
above a threshold value. Some conditions v/ill not permit 
location of the optimum if the step size is too small. The 
current step sizes, although generally too large, have per- 
formed well for the test craft above. The probable rationale 
is that the large step size permits a search over a wider, 
more extreme range of system values, but in the event that 
the step is too large, the step size will rapidly be 
reduced to an acceptable level. 

The lack of experimental data and analytical techniques 
has been noted in several areas. The nacelle design is a 
case in point. Although nacelles have been studied exten- 
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sively by the aeronautics industry, the application of 
nacelles to a water environment has been neglected and 
consequently, design considerations have not been defined. 

The design of cascaded corners has somewhat more experimental 
data for loss estimation, but the scatter of the data and 
the lack of definitive design information render elbow 
design a vague approximation at best. Three dimensional 
diffusers are likewise an unexplored domain. The effect of 
changes of ducting shape on losses is not treated in the 
literature for straight pipe. Loss estimation for changing 
shape in elbows is treated, but only for rectangular forms. 

It is apparent that the program results contain a great deal 
of information, but it is also evident that the interpretation 
of. the results can be improved by a more thorough study of 
runs for different craft, configurations, and missions. A 
study of the effects and interactions of such changes will 
provide a finer understanding of the import of the results. 
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CHAPTER 4 



CONCLUSIONS AND RECOMMENDATIONS 



A method of optimizing a waterjet propulsion system 
for hydrofoils has been developed and tested with a four 
hundred ton hydrofoil at a cruise speed of forty tv;o knots. 
Other craft have verified the importance of the selection 
of the gas turbine for a particular craft. It is 
recommended that a version of this program be adapted to 
determine the optimum displacement for a given engine. 

Although originally considered important, the amount 
of diffusion in the strut should be the minimum required to 
avoid cavitation in the ducting or the minimxm required by 
the pump, v/hichever is greater. The design of the nacelle 
need not be at the limits of cavitation, but should be 
integrated into the total system design to provide a lesser 
weight propulsion system. It is recommended that a study of 
nacelle design be initiated to determine the interaction 
of the nacelle and the propulsion system. 

The initial values of the system parameters may affect 
the operation of the optimization procedure. A determination 
of the cause of this lack of stability should be made. The 
magnitude of the step sizes of the parameters has little 
effect on the program provided they are sufficiently large. 

A general lack of data and design information is 
prevalent for fluid flow in closed conduits. This is 
especially noted for ram inlet design, cascaded corners, 

three dimensional diffusers, and changes of pipe shape. 
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The program provides a unique facility for the 
preliminary designer to raoidly ascertain the optimum 
water jet propulsion system design and performance. Many 
results should be studied to deduce the complete real 
world interpretation of a computer generated design. 
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APPENDIX A - OPTIMIZATION 



The selection of an optimization procedure represents 
a key consideration in the overall design method. For one 
thing it may require computations to be made solely for its 
usage, as in gradient methods. Then, too, it may require 
a supervisory program to tell it what to do, as in random 
search techniques, or an adaptive direction capability may 
be self-contained (ref. 18) . It is not the intent of this 
investigation to find or invent the "optimum" optimization 
procedure, but the significance of its usage should not be 
minimized. 

The general consideration of hov; to optimize has a 

wide variety of tools to aid in the achievement of the 

answer. One of the more common methods utilized is some 

variant of the gradient search (ref. 19) . The rationale 

is an intuitive one, based on the notion that if the goal 

is the attainment of a maximum, for example, the quickest 

way to go is in the direction of the steepest slope. This 

is not unlike a blind man's method of climbing a moxmtain. 

Like the blind man, the gradient search Ccin encounter 

difficulties by following this philosophy. The assumption 

of a finite slope in the area of interest not only requires 

the absence of steps on the hyper surf ace , but imposes the 

more stringent limitation that the slope be calculable and 

containable within the computer. This somev/hat limits the 

types of systems that can be optimized by this technique and 
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and exacts a stiff penalty for some numerical inaccuracies. 
Further, assuming that the system to be optimized permits 
slope calculations, it is not known a priori that such 
calculations provide sufficient information to compensate 
for the time spent in making them. In the present system 
design, discrete steps are virtually guaranteed due to the 
different pumps and gearboxes used. Since computation of 
the propulsion system weight is not done by means of an 
analytic function, slope calculations become as lengthy as 
the weight calculations themselves. For these reasons, 
gradient search techniques are rejected as a possible 
optimization tool for this design problem (ref. 20). 

On the opposite end of the optimization spectrum lies 
the random search. The principal advantages of a random 
search are the ease v;ith which it m.ay be understood, hence 
used and programmed, and the simplicity of its computation. 

No slopes or step sizes need to be calculated or estimated, 
merely the values of the function at the points of interest. 
This is also its major disadvantage. 

The random search also has no greater likelihood of 

finding a local optimum than it does of finding the global, 

and, hence desired optimum. Once the capability of 

"directionability is gained, local optima are just as likely 

to be found as the global optimum as is the actual optimum. 

This is true of all directed searches. The only way to 

avoid this unpleasant result is to use random starting points 
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or to use a systematic search over the entire hypersurface. 

A unimodal surface, of course, does not present this 
problem. 

Since the random search utilizes no other additional 
information than tlie function value, it is highly in- 
efficient with the result of greater time spent in the 
search. All optimization is done in comparison to the 
previous best. Thus the sinplicity of its procedure is 
gained only at the sacrifice of the large number of evaluations 
that must be made to assure a reasonably high certitude of 
attaining the optimum. This becomes especially significant 
when the time of calculation of a single function value is 
lengthy. In the design and performance estimation of a 
propulsion system, a lengthy calculation is likely to be the 
case, in part due to the iterative loops required (refs. 21, 

22 ) . 

The solution is clearly somewhere in the middle of 
these two approaches, utlizing the best features of each and 
minimizing their disadvantages. It should have the ability 
to discern direction without the requirement of additional 
calculations. One such method is the pattern search utilized 
in this design (ref. 23) . The basis of this method is akin 
to that of the gradient search but is modified in its ap- 
proach. Again based on intuition, the pattern search depends 
on the supposition that the search effort should be expended 

in the direction of improvement. Note that this is a 
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qualitative judgment rather than a quantitative one, as with 
gradient methods. The criterion for a sense of direction 
lies in the success of an evaluation, not in its measure 
of success, which in itself v;ould be somewhat arbitrary. 
Thus, while gaining adaptive information about the location 
of an optimum, the disadvantages of the gradient search are 
neatly bypassed and the random search is effectively uti- 
lized in a more efficient manner. 

Following the flow diagram of Figure 8, the pattern 
search uses the starting point as the initial basepoint. 

From the basepoint exploratory moves are made. That is, 
each coordinate is varied from the basepoint position the 
amovint of the step size and the function evaluated at that 
point. Each exploratory move is examined for possible 
improvement in the f motion value. If improvement cannot 
be made in any direction, the step size is decremented and 
the exploratory moves resumed from the basepoint. If a 
position is located with a better function value than the 
value at the basepoint, the second type of move, a pattern 
move, is made. A pattern move is a simultaneous change in 
all dimensions of length equal to the distance from the 
current basepoint to the present position. If improvement 
is made- in either the pattern move or the sequential explor- 
atory moves, another pattern move is made. Lack of improve- 
ment resets the basepoint at the position of the best 

fmction value and enters the cycle again with exploratory 
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moves from the current basepoint. 

This ability to seek, find, and establish a pattern 
gives this search an ability to move rapidly in the best 
direction. It also has the capability of shifting direc- 
tion or modifying the pattern when required. Step changes 
in the hypersurface present little problem since slopes are 
not utilized. With these features the pattern search can 
move quite readily along steep ridges or narrow valleys 
which are common "Waterloos" of other search methods, 
especially gradient methods. 

Perhaps it should be mentioned here that the intent 
is not to totally degrade the gradient search which has its 
own attributes and an important usefulness in its own right. 
It is not without shortcomings, however, that might present 
a serious detriment to the user who is looking for a 
generalized search method. 

The original form of the pattern search as proposed in 

reference 23 has been modified to conform to the needs of 

this design problem. An additional modification has been 

made that is considered an improvement on the initial 

scheme. During the exploratory move phase of the search, 

each dimension was initially incremented by a step size and 

the function evaluated. If the move was a success, the point 

and its associated value were saved and the next dimension 

explored. Lack of success meant decrementing the position 

by a step size and evaluating the function. Again a success- 
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ful move was followed by retaining the position and value and 
exploring the next dimension. An unsuccessful move was 
ignored and the next dimension incremented. The modification 
implemented is to vary the dimension first in the direction 
of the last improvement rather than simply incrementing it. 
Thus anticipating success, time would be saved by not 
checking a poor direction. The resulting advantage of this 
modification is, of course, very dependent on the function 
being evaluated, but is in keeping with the philosophy of 
the search and is accomplished at a very small cost of time. 
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APPENDIX B - NACELLE 



The details of this section may be found in reference 
9, from which the following comments are drawn. 

The selection of the system parameter values sets the 
ratio of the inlet diameter to the maximum external diameter, 
DIDM. This effectively limits that maximum amount of diffu- 
sion that may be carried out in the nacelle. The present 
design assumes that the maximum amount of diffusion allowed 
in the nacelle duct is utilized. The nacelle is first 
checked to ensure that it does not cavitate externally. 
Prediction of external cavitation rejects that selection of 
system parameter values. If cavitation is not predicted, 
the inlet is checked for cavitation at the point of maximum 
flow, which generally occurs in the take-off condition. If 
cavitation v/as predicted internally, auxiliary inlets are 
sized to reduce the inflow rate to just below the point of 
cavitation. The length of the nacelle is then suboptimized 
to be the length at which the total power loss due to drag 
and diffusion in the duct is minimized. The total head 
loss, made up of the lip loss and the diffuser losses may then 
be calculated. Likewise, the drag coefficient may be 
calculated from Hoemer's approximation for stramlined 
bodies (ref. 10), based on wetted surface, 

3 

Cd=Cf (l+1.5(^)^+7(I^m) 3) 
s Ln Lj^ 



- 51 - 



APPENDIX C - ELBOWS 



All elbows in the ducting are treated in identical 
fashion, although individually each elbow may be very 
different. The inlet and outlet areas for any bend are 
assumed to be the same. An elbow is considered fully defined 
by the knowledge of four quantities: width, depth, cross- 
sectional area, and radius ratio. The radius ratio, 
is the ratio of the radius of the centerline of the bend to 
the duct radius, as depicted in Figure 9. Since the flow 
rate, Q, has been previously determined, the first three 
items are sufficient to calculate the aspect ratio, shape, 
and the mean velocity of the fluid in the elbow. 

For Re ^ > 91, Ito's expression for the head loss 
coefficient in a circular bend is (ref. 11) ; 



K 



^=.OO241a0Re“® • ^’(| q) 



(C.l) 



a is a numerical factor deoendent on 0 and - . 

RO 

case of 0 = 90* and ^q<19.7 



For the 



a=0. 95+17.2 

RO 



(C.2) 



The loss coefficient then is 



Kt=0. 002410(0.95+17.2 Re"®* ^M|o) 

For a given duct and considering 6, V, and R constant, 
partial differentiation of with respect to RO leads to 
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5rO = 0 ^RO^opt " 

Cl 

Since K. = --^ + CoRO for these conditions, this result 
t RO 

verifies the intuitive feeling that a minimum exists due to 
a trade-off between friction and bend losses as the number 
of guide vanes in the elbov/ varies, as depicted in Figure 10. 

If the given radius ratio is greater than one, the 
losses in the elbow may be reduced by the use of splitters 
which aid in the prevention of separation. In this context, 
splitters are defined as guide vanes of unequal length, 
unequal spacing and possibly lanequal curvature (ref. 12). 
Turning vanes, on the other hand, are always equal in length 
and curvature, although the spacing between vanes is not 
necessarily equal. If the splitters are considered in order 
to sxabdivide the elbow into smaller elbows, the minimum 
loss will occur v/hen the ratio of the inside radius to the 
outside radius of any of the subdivided elbows is constant. 
Referring to Figure 6, this may be written as 



(ri) = 

A 

where n is the number of subdivided elbov;s, i.e. the number 
of splitters + 1. Figure 11 depicts a duct with one 
splitter. Clearly, 



r^ = R - RO 

r^ = R + RO 
-53- 



Hence, the desired radius ratio, 4.3, may be substituted for 



(roa~*~J^ia^ and the resulting equation solved for n, or. 



Ln 




2.12 Ln 




(C.4) 



The number of splitters determines the spacing. The head 
loss coefficient may be calculated for each subdivided 
elbow and the results averaged over the area which will 
yield the average head loss coefficient for the elbow. 

For R <1, equation C.4 has no meaning. For this case 



the total head loss coefficient was computed as the sum of 
the friction loss coefficient and the loss coefficient from 
Figure 12 which is an empirical graph for thin, circular arc 
turning vanes from reference 12. 

In either case the results are valid up to a Reynolds 
number below that expected in this design. The head loss 
coefficient including friction is then corrected for the 
higher Reynolds number by (ref. 13) 



RO 



K2 = ^ Ki 

Xi 
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where by curve fit, 

X = 1.0057-0. 16892Ln(Re x lO"®) +0 . 0145885 [Ln (Re x 10"®)]^ 

-0. 0016948 [Ln (Re x 10"®)]® (C.5) 
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APPENDIX D - STRUT 



The head losses in the strut diffuser are caused by 
friction eind the area expansion which thickens the boundary 
layer. The friction coefficient is computed in the same 
manner as for straight pipe. The expansion loss coefficient 
is 



C is dependent on the diffuser equivalent angle, 29, and 
is interpolated from the graph in Figure 13, which is for 
a two-dimensional, straight-walled diffuser. For the flov; 
area, the thickness to chord ratio of the strut is held 
constant while the area is increased. Hence the assumption 
of two dimensional diffusion is not strictly valid but is 
considered a conservative estimate. 

The strut dissipates energy externally by a drag term 
attributable to the body and a term resulting from the spray 
The drag coefficient from Hoerner (ref. 10) is based on the 
planform area and is approximated by 



The spray drag coefficient, likev/ise based on planform area 
is 



This approximation is valid for Froude number, based on the 
chord length, greater than three (ref. 14) . 




(D.l) 



Cdo = 2Cf^[l+2(t/c)+60(t/c) '*] 

s 



(D.2) 



Co = .03(t/c) 

'^SD 



(D.3) 



- 56 - 



APPENDIX E 



Pump Inlet Piping 

The ducting from the hull elbow to the pump inlet is 

a straight pipe, a 90® elbow and a transition piece. The 

length of the straight pipe depends on the number of prime 

movers. In equation form. 

Length = Beam 

Number of prime movers + 1 

The transition piece is shaped according to the number of 
pumps under consideration. Since there are two struts 
ducting v/ater into the hull, assuming a canard foil 
arrangement, and only one, two, or four pumps are consi- 
dered, the transition piece v;ill be a junction, straight 
pipe, or diverging branch. 

All straight pipe loss calculations are made using the 
standard Darcy-Weisbach equation 

h = fL (E.l) 

The Moody friction factor, f, is computed using the closed 
form approximation of Techo et al (ref. 15) , 

f = [. 86859 Ln (Re/1.964Ln(Re) - 3.8215))]“2 (E.2) 

R^ is based on the mean hydraulic diameter which is defined 
as 

Dg = 4A (E.3) 

P 

This approximation for the friction factor is accurate to 

within .1% over the range of Reynolds number from 10^ to 
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2.5 X 10^. 



For both the junction and divergence, the average 
velocity is assumed to be the same before and after the 
transition piece and the piping is symmetrical about the 
single pipe as shown in Figure 14. This requires that the 
smaller piping cross-sectional area is 1/2 that of the 
larger piping. 

The mixing loss coefficient in the junction, without 
friction loss, is confuted from 

Kj = 1 + X - 2cos (f (4>) ) (E. 4) 

where f((^) = 1.4(}) - .005834)^ 

and ({) is the angle the small branch makes with the larger 
duct (ref. 16). In the cases examined for this design, 4> 
v/as set to zero resulting in only friction losses in the 
junction. 

The divergence loss factor X was computed from Figure 
15 v/hich is adapted from reference 16. The angle of the 
divergence, of course, depends on the pump length. It 
should be noted that the divergence configuration requires 
another elbow, equal in angle to that of the divergence, 
to enter the punp. 

In all three transition pieces the size was set as 

required, depending on the pump length. If the pump's 

length requires a zero or negative transition piece 

length, the transition piece is neglected and the losses 

computed only for the athwartships length and the pximp elbow 
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which turns the ducting aft. 
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APPENDIX F 



Pump 

The details of this appendix may be found in reference 

17, from which the following comments are drawn. 

There are two basic types of purtps considered in this 

design; a multistage, double impeller centrifugal design 

and an axial pump with an inducer impeller. The variations 

occur in the number of stages used. The pump design 

utilizes the flow rate cind Thoma's criterion is defined as 

a = NPSH (F.l) 

Pump Head 

where the net positive suction head, 

NPSH = v2 + ha - h^ -helev " (F.2) 

2g 

and 

Pumo Head = v? - + h^i^„ + Ah^ (F.3) 

e- exev t 

The axial pump's first stage is an inducer impeller 
with up to two additional stages available if the amount 
of head is sufficiently high. Assumption of a specific 
speed yields the design pump speed and approximation 
to the off-design pump speed, which in this case is the 
cruise condition since the pump is basically designed 
for take-off conditions. The RPM of the prime mover deter- 
mines the gear ratio and calculation of the shaft horse- 
power permits an estimation of the gear weight. Finally, 

a fuel weight may be calculated for the range specified in 
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the cruise condition. 

The centrifugal pump is designed in the same v;ay except 
that the number of double suction impellers is increased, 
thus increasing the RPM, until the combined weights of the 
pump, dry and V7et, gearbox and fuel are at a minimum. 

The two resulting pump designs are then compared on 
the total pump, fuel, and gearbox weights and the least 
weight system is selected for the propulsion plant. 
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APPENDIX G 



Nozzle 

For the purpose of calculation the nozzle was consid- 
ered a straight-walled nozzle. Assuming the losses are of 
the .form 

h„ = f L (G.l) 

D 2^ 



and noting the friction factor is approximately constant 

and the velocity essentially uniform, this may be written 

Ah = f AL 
^ D 2g 

or 

L-; 2 

h_ = f 1 V dL 

0 D 2? 

Substituting D = + Ltcin(\j;) 

V = 4Q 

ttD^ 



results in 



h 



n 



16 

f Lj/2 [l-(i+Dt) •*] 
°t“^j TJJ 




(G.2) 



which is used for the head loss estimation in the nozzle 
when D > D- >1, (ref. 2). When D. <1, the head loss 

t j 3 

is assumed to be 1.5% of the dynamic head. 

In the event that the distance of the pump exit from 

the transom is greater than 2 times the height of the pump 

centerline above the keel, the nozzle is assumed to exit 

through the underside of the craft. The latter distance 
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is that if the nozzle were directed along an axis at 45® 
from the vertical. 

In every instance, however, the nozzle is assumed to 
be depressed at the optimum angle to reduce the lift 
required which would allow the use of smaller foils and 
hence reduce drag. As before, 

( 1 ) 

Rewritten this is also 

T = (W - T ) D (G.3) 

L 



T = pQ(VjCOS(6) - V^) 



where 

= pQVjSin(B) 

Equating equations (1) and (G.3), and solving for gives 

V. = (pQVq + WD) [cos (3) + Dsin (3)]“^ 

3 L L 

The minimum energy input to the jet occurs when is at a 
minimum, or, 

= 0 => (pQVq + WD)(Dcos(3) - sin (3)) 

33'^ L L = 0 

(cos (3) + Dsin (3) ) ^ 



=> 3 = TAN 



-1 



‘ E> 



(G.4) 

which is the optimum angle for nozzle depression (refs. 4,6). 
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RETURN TC CALLING 
PROGRAM 



LIST OF SYMBOLS 



Label 

FARMS 



DRAG 



FLOW 



USED IN COMMON 

Variable 

Name Subscript 



VJVO 

VIVO 

DIDM 

TDRAG(I) I 

1 

2 

3-5 

STRTD(I) 

POD(I) 

SPRAY (I) 

REST (I) 

VO (I) 

TRIM (I) 

Q(I) 
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Description/Usaqe 

System parameters 

Jet velocity ratio 

Inlet velocity ratio 

Inlet diameter to 
maximum nacelle diameter 
ratio 

Drag estimates 

Total craft drag at 
condition I, pound force 

Cruise condition 

Take-off condition 

Evaluation conditions 

Strut drag at condition 
I , pound force 

Nacelle drag at condition 
I , pound force 

Spray drag at condition 
1/ pound force 

TDRAG (I) -SPRAY (I) -STRTD (I) 
-POD(I)/ pound force 

Craft velocity at condition 
I, foot per second 

Craft trim angle at con- 
dition I, degree 

Flow characteristics 

Volume flow rate at con- 
dition I, cubic foot per 
second 



Label 



Variable 

Name 



Subscript Description/Usage 



BLEW 



CHARS 



AIN 

AJET 

AREA(ICOMP) 

VJ(I) 

VI (I) 



XK(IELB) lELB 

1 

2 

3 

4 

RO (lELB) 
THATA(IELB) 

WIDTH 



Inlet area, square foot 

Jet area, square foot 

Ducting cross-sectional 
area at station ICOMP, 
square foot 

Jet velocity at condition 
I, foot per second 

Inlet velocity at con- 
dition I, foot per second 

Elbov; data 

Radius ratio of elbow 
lELB 

Strut elbow 

Hull elbow 

Pump elbow 

Divergence elbow 

Duct radius at elbow 
lELB, foot 

Angle of bend of elbow 
lELB, degree 

Width of duct at elbow 
inlet 



DEPTH 



Depth of duct at elbow 
inlet 



TYPE (3, 1 TYPE) 



Contains name of elbow 
shape ITYPE 



WGTS (LS,LC) 
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System characteristics 

Weight of component LC, 
pound 



Label 



Vari able 
Name 



Subscript Description /Us age 



LS 



1 


Structure 


2 


Water (or fuel) 



LC 



1 


Nacelle 


2 


Strut elbow 


3 


Strut diffuser 


4 


Hull elbow 


5 


Athwart ships length 


6 


Pump elbov; 


7 


Transition piece 


8 


Pump 


9 


Nozzle 


10 


Reduction gear 


11 


Fuel 


12 


Prime mover 


13 


Lift from nozzle 
depression 


14 


Total weight 


15 


Spare location 


CGS (LG,LC) 


Centers of gravity of 
component LC, excluding 
fuel, gearbox, or prime 
mover, foot 



LG 



1 
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Structure vertical center 
gravity (from keel) 



1 



Label 



H20 



Variab le 

Name Subscript 

2 

3 

4 

DELH (I, I COMP) 

I COMP 
8 
9 

10-15 

CGSX 

CGSZ 

TEr^p 

PV 

RHOW 



Description /Usage 

Structure longitudinal 
center of gravity 
(from transom) 

Water vertical center of 
gravity, (from keel) 

VJater longitudinal center 
of gravity, (from transom) 

Total head loss up to and 
including component ICOMP 
at condition I, foot 

Same as defined except 

Nozzle head loss only 

Total system head loss 
(including elevation) , 
DELH(I,9)=DELH(I,8)+ 

DELH (1,7) 

Spare locations 

Total longitudinal center 
of gravity of system 
excluding gearbox, prime 
mover and fuel, from 
transom, foot 

Total vertical center of 
gravity of system exclud- 
ing gearbox, prime 
mover and fuel, from keel, 
foot 

Sea water (3.5% salinity) 
properties 

Temperature, degree 
Fahrenheit 

Vapor head, foot 

Density, pound force 
second squared per foot 
to the fourth 
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% 

'f 








pr" 

i . * 



d 






Label 



Variable 

Name 



TOLER 



PSUB 



Subscript 

GNU 

HA 

DELTA 

GERAT(NSTG) 

SHP (I,NSTG) 

RPM(I ,NSTG) 

PERF(L,IENGN) 

L 

1 

2 

3 

4 

5 

ETAP (I ,NSTG) 
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Description /Us age 

Viscosity, foot 
squared per second 

Atmospheric head, foot 



Check for optimum system 
Pump system 

Gear ratio required for 
pump NSTG 

Shaft horsepov/er re- 
quired at condition I for 
pump NSTG, horsepov;er 

Axial speed of pump NSTG 
at condition I, revo- 
lution per minute 

Prime mover lENGN 
Characteristics (L) 



Maximum normal horsepower 
at design speed, horse- 
power 

Maximum intermittent 
horsepower at design 
speed, horsepower 

Specific fuel consump- 
tion (SFC) at design speed 
and maximum normal horse- 
power, pound fuel per 
horsepower hour 

Design speed, RPM 

Prime mover weight, with- 
out auxiliaries, pound 

Efficienc^^ of punp NSTG 
at condition I 



Label 



Variable 

Name 



Subscript Description/Us age 



I 

7 

6,8 

SHIP 

DISP 

RAl'JGE 

BEAM 

HS 

HE 

HCL 

XLS 

XLPE 

XLP 

NACELLE 

DRAT 

DM 

AI 

AIAUX 

ELEXT 

ELENT 
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Same as defined except 
Reduction gear 
Spare locations 
Craft characteristics 
Craft displacement, pound 
Endurance, nautical mile 
Beam, foot 

Depth of sx±>mergence of 
foil, foot 

Height of pump centerline 
above keel, foot 

Height of pump centerline 
above mean water, foot 

Distance of centerline 
of strut root from tran- 
som, foot 

Distance of pump exit 
from transom, foot 

Length of pump, foot 

Nacelle characteristics 

Diameter ratio, DI/DM 

Maximum external diameter, 
foot 

Inlet area per nacelle, 
square foot 

Auxiliary inlet area per 
nacelle, square foot 

Length of forebody, foot 

Length of lip, foot 



Label 



CONST 



STRTC 



INDEX 



Variable 

Name Subscript 

ELAUX 

ELDIF 

ELN 

PI 

G 

RHOD 

TC 

T 

C 

T1 

Cl 

CFM 

lEVAL 

<0 

=0 

>0 



Description /Usage 

Length due to auxiliary 
inlet/ foot 

Length of diffuser, foot 
Length of nacelle, foot 
Constants 
3.14159265 

Acceleration of gravity, 

32. 174, foot per second 
squared 

Density of steel, 480, 
pound per cubic foot 

Strut characteristics 

Thickness to chord ratio 

Thickness at root, foot 

Chord at root, foot 

Thickness at tip, foot 

Chord at tip, foot 

Chord at flying water- 
line, foot 

Indices for program 
control 

Design/evaluation index 

Evaluate at |IEVAL| points, 
no design 

Design at cruise and take- 
off conditions only 

Design at cruise and take- 
off conditions and evaluate 
at |IEVAL| points 
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Label 



Vari able 

Name Subscript 

lEQPT 

=0 

>0 

ISTRT 

1 

3 

NUIIB 

2 

2+ I ievalI 

lENGN 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
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Des crip tion/Us age 

No equipments or config- 
uration entered; total 
system design/evaluation 

Specific equipments or 
configuration entered 
(not yet implemented into 
program) 

Initial program condition 
Design and/or evaluate 
Evaluate only 
Final program condition 
Design only 

Design and/or evaluate 
Prime mover type 
TF 35 
TF 40 



Proteus , 


1500 


rpm 


Proteus , 


1000 


rpm 



Tyne lA 
Tyne 1C 
FT12A 



LM 1500 
LM 2500 
FT4A-2C 
FT4A-12 



FT4C-2 



1 

% 

I 

I 



I 



Variable 
Label Name 



I TYPE 



I COMP 



NPUMP 

NGT 

ITABL 

CDRAG 

CSTRT(I) 
CPOD(I) 
CSPRY (I) 



Subscript 



1 

2 

3 

4 

1 

2 

3 

4 

5 

6 

7 

8 
9 



Description/Usage 

Elbow shape 

Ellipse 

Circle 

Rectangle 

Square 

Component index 

Nacelle outlet 

Strut elbow 

Strut diffuser exit 

Hull elbow 

Athv/artships length 

Pump elbow 

Transition piece 

Pump inlet 

Nozzle throat 

Number of pumps 

Number of gas turbines 

Interpolation parameters 

Computed drags 

Computed strut drag at 
condition I, poxind force 

Computed nacelle drag at 
condition 1 , poxind force 

Computed spray drag at 
condition I, pound force 



WEGT 
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Pump, gear and fuel weights 



Label 



Varicible 

Name 



Subscript Description/Us age 



PUMM 



XWD (NSTG) 



Pump NSTG dry weight, 
poxand 



xm<l (NSTG ) 
XWG(NSTG) 
XWF (NSTG) 

QQ(I) 

DIS (NSTG) 
D2S (NSTG) 



Pump NSTG wet weight, 
pound 

Gearbox weight associated 
with pump NSTG, pound 

Fuel weight associated 
with pump NSTG, pound 

Pump characteristics 

Flow rate per pump at 
condition I, cubic foot 
per second 

Inlet tip diameter of 
pump NSTG, foot 

Exit tip diameter of pump 
NSTG, foot 



XNS (NSTG) 



Specific speed of pump 
NSTG, cfs units 



SM(NSTG) 



Suction specific speed 
of pump NSTG, cfs units 



PLP (NSTG) 



Length of pump NSTG, foot 



NSTG 



1 

2 

3 

4 



Indicator of pvimp type 

Axial pump with single 
inducer impeller 

Axial pump with inducer 
and one impeller stage 

Axial pump with inducer 
and two impeller stages 

Centrifugal pump with max 
imum of ten parallel 
double suction impellers 
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5 



Spare location 




• "! 



' 














Variable 
Label Name 



SHI (NSTG) 



XIM 



HEAD 

HPP (I) 
HSV(I) 
THOM (I) 
PHI (NSTG) 
WF 
WG 

irkicic^ 

I FUEL 



Subscript Description /Us age 

Head coefficient of pump 
NSTG 

Number of parallel double 
suction impellers for 
centrifugal pump 

Number of impellers, not 
including inducer, for 
centrifugal pump 

Pump flow characteristics 

Pump head at condition I, 
foot 

Net positive suction head 
at condition I, foot 

Thoma's cavitation index 
at condition I 

Flow coefficient of pump 
NSTG 

Working variable for fuel 
weight, pound 

Working variable for gear 
weight, pound 

Un labeled common variables 

Logical variable for fuel 
calculation 

.TRUE. Make fuel calculation 

.FALSE. Do not make fuel calcula- 

tion 
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VARIABLES USED IN H20JT 
NOT IN COMMON 


Vari ab le 


Description 


FARM (3) 


VJorking system parameters 


ENGN(3,IENGN) 


Contains name of engine lENGN for output 


EHP (I) 


Effective horsepower at condition I, horsepowe 


DEL (3) 


Initial step sizes of system parameters 


DELMIN(3) 


Minimum step sizes of system parameters 


VK(I) 


Craft speed at condition I, knot 


IPRNT 


Print data set reference number 


DISPL 


Craft displacement, long ton 


IN 


Estimate of minimum number of gas turbines 
required to power craft 


MAX 


Estimate of maximum number of gas turbines 
needed to power craft 


XJ 


Working variable for number of gas turbines 


IK 


Working variable for craft condition I 


IL 


VJorking variable for IN 


KM 


Minimum number of pumps to be used for XJ 
gas turbines 


K1 


Maximum number of pumps to be used for XJ 
gas turbines 


WEIGT 


Propulsion system weight, long ton 


IM 


VJorking variable for prime mover character- 
istic L 


DISPL 


Craft displacement, long ton 


I 


Index for do loops 


IJ 


V'Jorking variable for check on power adequacy 


J 


Index for do loop 
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Variable 



Description 



I SAVE 
N 

XfllN 

MK 



Working variable for MAX 

Working variable for output of prime mover 
type 

Factor for minimum acceptable prime mover 
SHP operation 

Working variable for prime mover character- 
istic L 
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LIST OF VARIABLES USED IN FCT 
NOT IN COMMON 



Variable 


Description 


APOD (I) 


Difference betv;een present and previous 
nacelle drag calculations at condition I, 
pound force 


ASTRT(I) 


Difference between present and previous strut 
drag calculations at condition I, pound 
force 


ASPRY (I) 


Difference between present and previous spray 
drag calculations at condition I, pound 
force 


C(I) 


Margin factor for thrust at condition I 


TSUM 


Working variable for propulsion system 
weight, pound 


DSIG 


Working variable for SIGIIA 


TOTAL 


Total head less vapor head, foot 


SIGMA 


Head above vapor head, foot 


SUMZ 


Working variable for moment in vertical plane 
due to weights, foot pound force 


SUI>IX 


Working varicible for moment in longitudinal 
plane due to weights, foot poiand force 


SUM 


Working variable for system weights, pound 


ZSUM 


Total moment in vertical plane due to weights 
foot pound force 


XSUM 


Total moment in longitudinal plane due to 
weights, foot pound force 


SPD 


Working variable for jet velocity calcula- 
tion, foot per second 


J 


Working varicible for NUMB 


JDRAG 


Indicator of acceptable drag accuracy 


1 


New drag calculation is within 5% of previous 
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1 



Variable Description 





calculation 


2 


New drag calculation is greater than 5% of 
previous calculation, redesign strut and 
nacelle 


ANGLE 


Optimum nozzle depression angle, radian 


COEF 


Factor for nacelle drag calculation 


0 


Strut must be resized due to cavitation 


1 


Use nacelle design 


KOUNT 


Working variable for craft condition I 


I 


Craft condition 


FARM (K) 

K 


Working variable for system parameters 


1 


VJVO 


2 


VIVO 


DDRAG(I) 


Previous acceptable total drag calculation at 
condition I, pound force 


DSPRY(I) 


Previous acceptable spray drag calculation at 
condition I, pound force 


DSTRT(I) 


Previous acceptable strut drag calculation at 
condition I, pound force 


DPOD(I) 


Previous acceptable nacelle drag calculation 
at condition I, pound force 



LIST OF VARIABLES USED IN NAGEL 
NOT IN COMMON 



Variable 


Description 


ZK 


Percentage of auxiliary inlet area permitting 
flow 


SPO 


Static pressure at inlet, pound per square 
foot 


PVP 


Vapor pressure, pound per square foot 


SIGTV 


Incipient turning vane cavitation number, 
referenced to diffuser exit velocity 


JNUMB 


Working variable for NUI^tB 


I 


Craft condition, also index for do loops 


CPEX 


Peak external pressure coefficient 


K 


Working index for do loop for ELEXT/DM 
values 


J 


Working index for do loop for VIVO values 


XD 


Length of forebody to maximum diameter ratio, 
ELEXT/DM 


DIDMX 


Maximum permissible diameter ratio 


QI 


Flow rate per strut, cubic foot per second 


DI 


Nacelle inlet diameter, foot 


CP IN 


Peak internal pressure coefficient 


QIN 


Take-off flow rate per nacelle through inlet, 
cubic foot per second 


QC 


Take-off flow rate per nacelle, total, cubic 
foot per second 


QAUX 


Auxiliary flow rate required to avoid 
cavitation, cubic foot per second 


KDEX 


Counter for iterations on diffuser length 


VR 


Velocity ratio at take-off condition, based 
on flow through inlet 
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Variable 


Description 


VI 2 


Inlet velocity at take-off condition, 
based on flow through inlet 


PRL2 


Pressure recovery coefficient of lip at 
take-off 


PTI 


Inlet stagnation pressure immediately aft 
of lip at take-off, pound per square foot 


SPI 


Inlet static pressure immediately aft of lip 
at take-off, pound per square foot 


VIAUX 


Inlet velocity in auxiliary inlet, foot per 
second 


PRAUX 


Pressure recovery of auxiliary inlet 


PTAUX 


Stagnation pressure inside auxiliary inlet, 
povind per square foot 


DYP 


Ilet dynamic pressure immediately aft of lip, 
pound per square foot 


PC 


Average inlet stagnation pressure of combined 
flow, pound per square foot 


QDIF 


Working variable for diffuser flow rate, 
cubic foot per second 


PHI 


Equivalent angle of forebody, radiain 


PHS 


Sine of equivalent angle of forebody 


X 


Length of auxiliary inlet, foot 


D2 


Diffuser exit diameter, foot 


D1 


Inlet diameter, foot 


ELMAX 


Maximum permissible length of diffuser, foot 


ELMIN 


Minimum permissible length of diffuser, foot 


II 


Working variable for craft condition I 


EL 


Working variable for diffuser length, foot 


DEL 


Working variable for change in diffuser length, 
foot -89- 
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Vari ab le 


Description 


ELD 


Working variable for nacelle length, foot 


ELL 


Working variable for nacelle length, foot 


ELF AC 


Working variable for excess in nacelle length 
due to diffuser, foot 


DDM 


Average diffuser diameter, foot 


XKT 


Form loss coefficient of diffuser 


REL 


Reynolds number, based on nacelle length 


RED 


Reynolds number, based on inlet diameter 


DL 


Ratio of maximum external diameter to nacelle 
length 


CDRG 


Corrtputed drag coefficient 


AJJGL 


Equivalent half angle of diffuser, degree 


CDIF 


Diffuser expansion factor 


ROW 


Power loss due to drag and duct loss of 
diffuser, horsepower 


POWI 


Previous power loss calculation for diffuser 
length, horsepower 


EM 


Factor in wetted surface calculation 


AEXN 


Wetted surface area, square foot 


REND 


Reynolds number, based on inlet diameter 


DDIF 


Total pressure loss in diffuser, pound per 
square foot 


FLOSS 


Total pressure loss in nacelle, pound per 
square foot 


VAOUT 


Average exit velocity, foot per second 


SQUAR 


Factor in critical velocity calculation 


VCRIT 


Critical velocity in strut elbow, foot per 
second 
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Variable 


Description 


VMAX 


Maximum velocity at nacelle exit, foot per 
second 


RENL 


Reynolds number, based on nacelle length 


NL(2) 


Interpolation parameters 


ML (2) 


Interpolation parameters 


KL(2) 


Interpolation parameters 


JL(2) 


Interpolation parameters 


IL(2) 


Interpolation parameters 


VRT(6) 


Data array of velocity ratios 


XDT(IO) 


Tabulated forebody length to inlet diameter 
ratios 


PRLT (6) 


Tabulated lip pressure recovery coefficients 


SIGI (I) 


Free stream cavitation index at condition I 


PTO(I) 


Stagnation pressure at craft condition I, 
pound per square foot 


XDTT(J) 


Tabulated ELEXT/DM ratios for trim angle J 


CDUMX (K) 


Dummy array of peak pressure coefficients for 
VIVO K 


CDUMY(L) 


Dummy array of velocity ratios for ELEXT/DM 
L 


CD(I) 


Drag coefficient at craft condition I 


DIDMT(J) 


Tabulated DI/DM values for ELEXT/DM J 


VRTT(J) 


Dummy array of velocity ratios for angle of 
attack (J) , internal 


DLIP (I) 


Lip loss coefficient at craft condition I 


QO(I) 

VELR(I) 


Free stream dynamic pressure at craft condition 
I, pound per square foot 

Inlet velocitv ratio at craft condition I 
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VELR(I) 



Variable 


Description 


VRTEX (I) 


Same as VRTT(I) but external 


VRMAX(I) 


Maximum permissible velocity ratio at craft 
condition I 


FRICT 


Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 
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VARIABLES USED IN ELBOW 
NOT IN COmON 



Variable 
SHAPE(3,ITYPE) 
THETA (J) 
XLOSS(J) 

ROA(IO) 

RE 

FRICT 

CORR 

RExMAX 

KOUNT 

lELB 

1 

2 

3 

4 

FACTR 

RIN 

ROUT 

RATIO 

XN 



Description 



Contains name of shape ITYPE 
Data array of elbow angles 

Data array of elbow loss coefficients with 
thin, circular arc turning vanes 

Outside radius of splitters, foot 

Reynolds number of duct 

Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 

Function statement ,• calculates head loss 
correction factor for different Reynolds 
numbers 

Maximum Reynolds number permitted for 
splitter loss equation 

Working variable for craft condition I 

Index indicating which elbow is being 
designed/evaluated 

Strut 

Hull 

Pump 

Divergence 

Factor used in splitter loss calculation 

Inside radius of bend, foot 

Outside radius of bend, foot 

Desired radius ratio, 4.3 

Number of subdivided elbows required to 
achieve RATIO=4.3 
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Variable 


DescriDtion 


N 


Nxmber of subdivided elbows used 


N1 


Number of splitters corresponding to N 


SUM 


Working variable of svm of head times 
subdivided elbow area, foot cubed 


RIA 


Inside radius of subdivided elbow, foot 


V 


Average velocity, foot per second 


HGT 


Height of subdivided elbow, foot 


AA 


Equivalent cross-sectional area of subdivided 
elbow based on HGT, square foot 


RAD 


Equivalent radius of subdivided elbow, foot 


XCORR 


Ratio of head loss correction factors for 
differing Reynolds numbers 


XKT 


Head loss coefficient of subdivided elbow 


DIAM 


Equivalent diameter, foot 


VOLV 


Volume of splitters, cubic foot 


VOL 


Volume of splitters and elbow structure, 
cubic foot 


AREAl 


Duct area per elbow, square foot 


IJ 


VJorking variable for shape determination 


IK 


Working variable for shape determination 


I 


Working variable for number of subdivided 
elbows N 


RATEO 


Ratio of inside radius to outside radius 
of subdivided elbows 


HEAD 


Average head loss in elbow, foot 
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LIST OF VARIABLES USED IN STRUT 
NOT IN COMMON 



Variable 


Description 


THET2 


Array of data of equivalent angle of diffuser 
degree 


EXPAN 


Array of data of expansion coefficient 


ARATO 


Area ratio of strut diffuser 


VL 


Local maximum velocity, external, foot per 
second 


SIGMA 


Local cavitation number, external 


ai 


Mean chord, foot 


DEIN 


Strut inlet equivalent diameter, foot 


WIDE 


Width of duct at strut exit, foot 


DEOUT 


Strut exit equivalent diameter, foot 


DEAVE 


Average equivalent diameter, foot 


STRT 


Vertical strut length, foot 


XLONG 


Actual strut length, foot 


STAl^ 


Arctangent of equivalent diffuser angle 


THETA 


Equivalent angle of diffuser, 20, degree 


ECOEF 


Diffuser expansion factor 


FORI-IL 


Diffuser expansion loss coefficient 


VELIN 


Average inlet velocity, foot per second 


VLOUT 


Average exit velocity, foot per second 


RES 


Strut Reynolds number, based on mean chord 


CDS 


Strut drag coefficient 


CDSP 


Spray drag coefficient 


P£ 


Duct Reynolds number, based on inlet velocity 
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Variable 



Description 



PIPED 


Duct friction loss coefficient 


TOTAL 


Total loss coefficient 


HEAD 


Head loss of diffuser, foot 


KOUNT 


Working variable for craft condition I 


HGT 


Elevation of strut, foot 


CGWS 


Vertical center of gravity of duct, foot 
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LIST OF VARIABLES USED IN JUNCT 
NOT IN COMMON 



Variable 


Description 


ALPHA 


Assumed junction angle of 0®, degree 


DIAM 


Duct diameter, foot 


FRCTL 


Friction loss coefficient for athwartships 
length 


KOUNT 


Working variable for craft condition I 


AJCT 


Length of fore and aft ducting to pump inlet, 
foot 


FRCTJ 


Friction loss coefficient for junction 


XLAMD 


Working variable for AMIXL calculation 


FLONG 


Length of athwartships ducting, foot 


AJCTL 


Total loss coefficient of junction 


AMIXL 


Mixing loss coefficient of junction 


V 


Average velocity, foot per second 


XPUMP 


Number of gas turbines + 1 


RE 


Reynolds number of duct 


BETA(J) 


Data array of junction angles, degree 


ALAMD(J) 


Data array of mixing loss coefficient 
corresponding to BETA(J) 


FRICT 


Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 
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LIST OF VARIABLES USED IN PIPE 
NOT IN Ca-IMON 


Variable 


De script! on 


KOUNT 


Working variable for craft condition I 


V 


Average velocity, foot per second 


DIAM 


Duct diameter, foot 


RE 


Reynolds number of duct 


APIPE 


Length of fore and aft pipe to pump inlet, 
foot 


FRCTL 


Friction loss coefficient for athwartships 
length 


XPUMP 


Number of gas turbines 


XKT 


Friction loss coefficient of fore and aft 
length 


XLONG 


Length of athwartships pipe, foot 


FRICT 


Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 
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LIST OF VARIABLES USED IN DIVRG 
NOT IN COMt^ON 



Variable 


Description 


FRICT 


Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 


KOUNT 


V7orking variable for craft condition I 


RE 


Reynolds number of duct 


DIVL 


Divergence loss coefficient, not including 
friction 


DWGTl 


Duct weight of divergence angle, pound 


DWGT2 


Duct weight of divergence length, pound 


DV'?GT3 


Duct weight of pump inlet angle, pound 


W'JGTl 


Water weight of divergence angle, pound 


t\TWGT2 


Water weight of divergence length, pound 


miGT2 


V7ater weight of pump inlet angle, poiind 


CGWXl 


Longitudinal center of gravity of divergence 
angle, from transom, foot 


CGWX2 


Longitudinal center of gravity of divergence 
length, from transom, foot 


CGWX3 


Longitudinal center of gravity of pump inlet 
angle, from transom, foot 


DIAM 


Duct diameter, foot 


FRCTL 


Friction loss coefficient for athwartships 
length 


DIVLC 


Total divergence loss coefficient 


XPUMP 


Number of gas turbines + 1 


ANGLE 


Divergence angle, also pump inlet angle, radian 


HEADL 


Total divergence head loss, foot 


FLONG 


Athwartships length, foot 
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Varicible 

XLONG 

V 

ADIV 
THETA (J) 
COEF(J) 



Description 

Fore and aft length, foot 

Average velocity, foot per second 

Divergence length, ADIV = X LONG/COS (ANGLE) , foot 

Data array of divergence angles, degree 

Data array of divergence loss coefficients 
without friction, corresponding to THETA (J) 
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LIST OF VARIABLES USED IN NOZZL 
NOT IN COMMON 



Variable 


Description 


KOUNT 


Working variable for craft condition I 


XLPS 


Dummy nozzle length, foot 


XFAC 


Check on whether nozzle exits through bottom 
or stem 


XPUMP 


Number of pumps 


Ai'IOZ 


Optimum nozzle depression angle, radian 


XLNOZ 


Nozzle length, foot 


DT 


Nozzle throat diameter, foot 


DJ 


Nozzle jet diameter, foot 


XCORR 


Nozzle head loss factor 


RE 


Reynolds number, based on average diameter 
and velocity 


AREAl 


Throat pipe area, square foot 


AJETl 


Jet pipe area, square foot 


QQ 


Flow rate per nozzle, cubic foot per second 


FRICT 


Function statement, calculates Moody friction 
factor for smooth pipe, based on Reynolds 
number 
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LIST OF VARIABLES USED IN PUMP 
NOT IN COiMMON 



Variable 


Description 


K 


Cruise condition indicator 


J 


Take-off condition indicator 


I 


Index for do loops 


XPUMP 


Number of pumps 


THOMI 


Lower limit of Thoma's criterion for single 
inducer axial pump 


DRAT 


Hub to tip diameter ratio 


XX LP 


Factor for pump length 


CW 


Weight coefficient 


QX 


Ratio of cruise flow rate to take-off flow 
rate 


CA 


Factor for RPM calculation 


CB 


Factor for RPM calculation 


HX 


Patio of design pump head to off design head 


RX 


Ratio of design to off design pump RPM 


PHRAT 


Ratio of off design to design flow 
coe f ficients 


XNGT 


Number of gas turbines 


ETAPP 


Product of pump and gearbox efficiencies 


HP 


Inducer head for one and two stage axial 
pump designs, foot 


THOMS 


Thoma's cavitation criterion for inlet to 
axial stage 


HHP 


Axial stage pump head, foot 


BETA2 


Exit blade angle, radian 
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Vari ab le 


Description 


XNNS 


Non-dimensional specific speed 


BD 


Impeller exit v;idth ratio 


CC 


Factor in flow coefficient calculation 


AA 


Factor in flow coefficient calculation 


IMPL 


Maximum number of impellers permitted for 
centrifugal pump, = 10 


M 


Working variable for nximber of centrifugal 
pump impellers 


N 


Working variable for NSTG 


PC(K, J) 


Inducer head curve coefficients 


PCA(K, J) 


Inducer plus axial stage head curve coefficients 


PCC(K, J) 


Centrifugal pump head curve coefficients 


XRPM(M) 


Working variable for off design RPM of 
centrifugal pump with M impellers, RPM 


XDl (M) 


Working variable for inlet tip diameter, DlS, 
foot 


XD2 (M) 


Working variable for exit tip diameter, D2S, 
foot 


RPK(M) 


Working variable for design RPM of centrifugal 
pump with M impellers, square foot 


XPUP (M) 


Working variable for centrifugal pump area 
with M impellers, foot squared 


YLP (M) 


Working variable for centrifugal pump length 
with M impellers, foot 


XERAT(M) 


Working variable for gear ratio required for 
centrifugal pump with M impellers 


APUP(NSTG) 


Inlet area of pump NSTG, square foot 


WRAT (NSTG) 


Weight ratio of pump NSTG, including pump dry 
and wet weight, gearbox and fuel 
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Variable 


Description 






WD(M) 


VJorking variable for XWD(NSTG) for 
centrifugal pump with M impellers 




VW(M) 


Working variable for XVTW(NSTG) 
fugal pump with M impellers 


for 


centri 


\WG<M) 


Working variable for XWG(NSTG) 
fugal pump v/ith M impellers 


for 


centri 
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LIST OF VARIABLES USED IN FUEL 
NOT IN C0M40N 



Variables 


Description 


CFS 


Constant for SFC calculation 


CA 


Cruise condition drag to lift ratio 


CD 


1 + total system head loss coefficient, 
based on jet velocity 


N 


Number of intervals endurance is divided into 
+ 1 


XN 


Number of intervals endurance is divided into 


TI 


Time to cover one range interval at constant 
VO(l), hour 


UK 


Index for shifting SFC curves 


XJ 


Factor to convert SFC if SHP is less than 70% 
of design SHP 


I 


Working index for N 


H 


Head required, foot 


SHPP 


Total thrust required, horsepower 


SHNG 


Total thrust required per engine, horsepower 


M 


Working index for N 


WT(ITIflE) 


Weight of fuel used in time increment ITIME, 
pound 


DIS (ITIME) 


Displacement at time increment ITIME, pound 


VJJ(ITIME) 


Jet velocity at time increment ITIME, foot 
per second 



- 105 - 



I 



#■ 




<3 



. l-ssi-S/'. ■ ’’ 




LIST OF VARIABLES USED IN PTTRN 
NOT IN COMNON 



Vari able 


Description 


PSI (N) 


Current basepoint coordinate of parameter N 


THETA (N) 


Previous basepoint coordinate of parameter N 


PHI (N) 


Present exploratory point coordinate of 
parameter N 


DEL(N) 


Current step size of parameter N 


DELMIN(N) 


Minimum step size of parameter N 


DIR(N) 


Last successful direction of parameter N 


SAVE (N) 


Working variable for PHI(N) 


S 


Working variable for function value 


SPHI 


Working variable for function value at PHI 
coordinates 


SPSI 


Current best function value at PSI coor- 
dinates 


NUMB 


Counter for minimum step size check 


RHO 


Step size change factor 


I CALL 


Indicator of current point move 


K 


Index for do loop 


I 


Index for do loop 


N 


Number of parameters of search 


SIGN (N) 


Directed step size of parameter N 
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LIST OF VARIABLES USED IN OUTPUT 





NOT IN COMMON 


Vari ab le 


Description 


IPRI'IT 


Print data set reference number 


J 


Working variable for condition I 


SUM 


Total duct head loss, excluding elevation, 
foot 


I 


Index for do loop 


DRATO 


Strut diffuser area ratio 


IK 


Index for number of elbows in system 


XNGT 


Number of gas turbines 


K 


Index for do loop 


L 


Index for implied do loop in output statement 


LL 


Working variable to point to correct head 
loss for output 


KLL 


Working variable to point to correct format 
statement for head loss output 


M 


Index for implied do loop in output statement 


NIMP 


Number of impellers in p\amp 


TFM 


Strut thickness at flying waterline, foot 


WTRAT 


Total propulsion system weight ratio 


WRATF 


Propulsion system weight ratio, excluding fuel 


ENGN(IENGN) 


Contains name of engine lENGN 


VJRAT (I) 


Jet velocity ratio at craft condition I 


VI RAT (I) 


Inlet velocity ratio at craft condition I 


HEADL(I,ICOMP) 


Head loss of component ICOMP at craft con- 
dition I 


CONDS(2,I) 


Label for craft condition I 
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Variable 



Description 



ELBWS ( 2 
PC(I) 
LABEL (5 



IK) Label for elbow IK 

Propulsive coefficient at craft condition I 
M) Labels for output 
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